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Chapter 1 


Free Oscillations of a Linear 
Oscillator 


1.2 Review of the Principal Formulas 


The differential equation of a free linear torsion oscillator: 
G+ 2p + wey = 0. (1.1) 


The frequency and the period of free oscillations without friction (at y << wo): 


D 2 
ge aes” (1.2) 
J Wo 


An oscillatory solution (valid at y < wo): 
y(t) = Ape cos(w it + ôo), (1.3) 


where the constants Ag and ôo are determined by the initial conditions (0), (0). 
The frequency wı of damped oscillations 


wi = ywi = 7. (1.4) 
An equivalent form of the general solution: 
y(t) = e7% (C coswit + Ssinwyt), (1.5) 


where the constants C and S are determined by the initial conditions. They are 
related to Ap and 6p: 


Ap = VC2 +82, — tandy = —S/C. (1.6) 
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In the case of weak damping (y < wo) 
w1 S wo — 77 /(2wo). (1.7) 
The decay time (during which the amplitude is reduced by the factor e ~ 2.72): 
T =y: (1.8) 
A non-oscillatory motion at y = wo: 
p(t) = (Cit + Co)e™™. (1.9) 


The quality factor Q of an oscillator: 


The number of oscillations, during which the amplitude is halved: 


ln2 Q 
Ni; = — Q = 0.22 Q = —. 1.11 
Dam Q =0.22 Q 153 (1.11) 
The total mechanical energy of the oscillator consists of elastic potential energy 
of the strained spring and kinetic energy of the flywheel: 


1 


E= Epot + Exin = 2 


1 
Dy? + z7% (1.12) 
The values of the potential energy and kinetic energy of the oscillator, averaged 
over a cycle, equal one another, each of them constituting one half the total energy: 


1 1 1 
(Epot) = (Ekin) = zË = P4 = 172040 (1.13) 


1.3 Questions and Problems with Answers and 
Solutions 


1.3.1 Free Undamped Oscillations 


1.3.1.1 The Initial Conditions and the Shape of the Plots. In the absence 
of friction a linear oscillator executes simple harmonic motion, which is charac- 
terized by purely sinusoidal time dependence of the angular displacement and of 
the angular velocity. 

(a) What initial conditions give rise to oscillations of cosine time dependence, 
of sine time dependence? Suppose that you want to get oscillations with the angu- 
lar amplitude of 90°. What initial angular displacement (0) = yo at zero initial 
angular velocity ¢(0) = 0 ensures the desired amplitude? 
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Free oscillations with a cosine time dependence of y(t) = po cos wot 
occur if we set some arbitrary nonzero initial displacement y(0) = 
(vo, and the initial velocity zero: (0) = 0. To ensure an amplitude 
of 90° in this case, the initial deflection pọ must be +90° (47/2 
rad). On the other hand, if we set a (nonzero) initial angular velocity 
(Y(0) = Q) and the initial displacement zero (y(0) = 0), we obtain 
subsequent natural oscillations that are described by a sine function: 
p(t) = (Q/wo) sin wot. 


(b) What initial angular velocity (0) = Q ought you to impart to the oscil- 
lator, at rest in the equilibrium position, in order to obtain the same amplitude of 
90°? Remember, that the initial angular velocity Q must be expressed for input in 
units of the natural frequency wo. Verify your answer with a computer experiment, 
using the appropriate initial conditions. 


If the oscillator is excited by an initial velocity (g(0) = Q) at the 
initial displacement zero (y(0) = 0), the amplitude of resulting os- 
cillations equals |Q|/wo. To ensure oscillations whose amplitude is 
90° (7/2 rad), we need to input Q = +(7/2)wo (£1.57 in units of 
the natural frequency wo). 


1.3.1.2 Maximal Deflection and Conservation of Energy. Imagine exciting 
an oscillator initially at rest in the equilibrium position by a push which produces 
an initial angular velocity Q = 2wo. 

(a) Calculate the angle Ymax of maximal deflection using the law of the con- 
servation of energy. 


To obtain this result from the law of energy conservation, we should 
set equal the initial kinetic energy Ekin = iJ Q? to the potential 
energy Epot = t DYZ ax at the maximal deflection Ymax. Thus we 
find the same value for the maximal deflection from the equilibrium 
position: Ymax = QV J/D = N /wo = 114.6°. 


(b) Verify your result experimentally. Note that the simulation program per- 
forms the numerical integration of the differential equation independently of con- 
servation laws, such as the conservation of energy. That is, these laws are not 
used in the program. If we set the initial angular velocity Q = 2wọ, and the initial 
displacement zero, then the amplitude (and the maximal deflection Ymax from the 
equilibrium position) equals Q/wo = 2 radians, or 360° /r = 114.6°. 


1.3.1.3 The Phase Trajectory and the Initial Conditions. Compare the 
motion of the representative point along the phase trajectory of a conservative os- 
cillator with the time-dependent plots of the angle of deflection and of the angular 
velocity. 

(a) How is the phase trajectory changed if you change the initial conditions? 


https://gioumeh.com/product/simulations-—of-oscillatory-systems-solution/ 


8 CHAPTER 1. FREE OSCILLATIONS OF A LINEAR OSCILLATOR 


When we change the initial conditions, we change the initial mechan- 
ical state of the system. In the phase plane, changing the initial con- 
ditions means changing the point from which the phase trajectory 
originates. Since for the system with given properties only one phase 
trajectory passes through a given point of the phase plane, the further 
motion of the representative point is defined uniquely. 


(b) Does the direction of the motion of the representative point along the phase 
trajectory depend on the initial conditions? 


The representative point always moves clockwise along a closed phase 
trajectory of a conservative linear oscillator, independently of the ini- 
tial conditions. 


(c) Is it possible that phase trajectories for different initial conditions coincide? 
If so, formulate the requirements for the coincidence. 


For a conservative system, phase trajectories corresponding to differ- 
ent initial conditions coincide if these initial conditions correspond 
to the same total energy. However, the representative point starts its 
motion at different points of this phase trajectory depending on the 
initial conditions. 


1.3.1.4 Elliptical and Circular Shape of the Phase Trajectory. 

(a) Prove analytically that the phase trajectory of a conservative linear oscil- 
lator is an ellipse with its center at the origin of the phase plane. What are the 
semiaxes of the ellipse? 


To obtain the explicit equation of the phase trajectory, i.e., the depen- 
dence of ġ on ọ for a given motion, we can combine the expressions 
for y(t) and ġ(t) in the general solution of the differential equation 
of motion in order to eliminate time t. Since sin? a + cos? a = 1, we 
obtain: 


=i (1.14) 


Thus the phase trajectory of a conservative linear oscillator is an el- 
lipse whose center is located at the origin of the phase plane and 
whose semiaxes are Ap and Aowo respectively. 


(b) Show that the elliptical shape of the phase diagram of a conservative linear 
oscillator follows immediately from the law of the conservation of the energy. 


For a conservative system, the relation between ¢~ and y that gives 
the equation of the phase trajectory, follows directly from the law of 
energy conservation. This means that actually we don’t need to solve 
the differential equation of motion in order to get the equation of the 
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phase trajectory. Indeed, equating the sum of kinetic and potential en- 
ergies at an arbitrary time to the constant value of total energy (which 
equals the potential energy at maximal deflection Ag), we find: 

1 1 

D? + Je? = =DAG. 

i ae ae a 

Since J/D = 1/w%, this is the equation of the same ellipse, Eq. (1.14), 
which has been obtained above from the general solution of the dif- 
ferential equation of motion of the system. 


(c) What scale on the axis of the ordinate (the angular velocity axis) of the 
phase plane produces a circular phase trajectory? 


The phase trajectory of a conservative linear oscillator is circular if we 
plot the angular displacement ọ along the abscissa axis of the phase 
plane in radians and the angular velocity ¢ along the ordinate in units 
of the natural frequency wọ. In other words, the phase trajectory is a 
circle if we plot along the ordinate axis the quantity ~/wo in radians 
using the same scale both for the abscissa and ordinate axes. In this 
case the semiaxes of the ellipse defined by Eq. (1.14) are equal. 


(d) Does the time interval during which the representative point passes along 
one loop of the phase trajectory depend on the initial conditions? 


The radius vector of the representative point in the phase plane rotates 
clockwise with the mean angular velocity wo which equals the natural 
frequency of the oscillator. One revolution of the radius vector is 
accomplished during one complete cycle of oscillations. Rotation of 
the radius vector is uniform if the phase trajectory is a circle. If the 
phase trajectory is an ellipse, the radius vector of the representative 
point rotates nonuniformly. However, the average angular velocity is 
always wo. 


1.3.1.5 The Phase Diagram and Energy Transformations. Compare the 
phase trajectory with the plot of potential energy versus the angle of deflection. 
The positioning of plots on the display screen (if you open the window “Phase di- 
agram”) is convenient for such comparison. Pay special attention to the positions 
of the extreme points (turnarounds) on the phase trajectory and in the parabolic 
potential well. For the initial conditions y(0) = yo, (0) = Q, what are the 
values of the potential energy and the kinetic energy at the extreme points and at 
the equilibrium position? 

What are the extreme deflection Ymax and the maximal angular velocity wmax 
of the flywheel? 


For the initial conditions y(0) = yo, Y(0) = Q, the total energy is 


1 1 
Et = —Dy? + -JQ . 
tot 5) Pota 
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At the extreme deflection, Ymax, the kinetic energy is zero, and the 
total energy equals the potential energy of the strained spring: 


1 1 1 
z PPinax _ zP% + zI hence Ymax = 1/ 9% + ?/uz. 


At the equilibrium position y = 0, potential energy of the spring 
becomes zero, and all the energy is kinetic energy of the flywheel, 
rotating with the maximal angular velocity Wax: 


1 1 1 
z Tmax = Ate z PPO hence wmax = 4/2? + weye. 


1.3.1.6 The Shape and the Frequency of Energy Oscillations. Consider the 
plots of the time dependence of kinetic energy and potential energy. 

(a) What can you say about their maximal and average values? Compare these 
plots with the plots of the angular displacement and the angular velocity. 

(b) At what frequency do the oscillations of each kind of energy occur? What 
are the limits (the extreme values) and the mean (averaged over a period) values 
of each kind of energy in these oscillations? 


(a, b) At free oscillations in any conservative system the maximal 
values of each kinetic and potential energy coincide with the total en- 
ergy Erot, whose value remains constant during the oscillations. In 
a linear conservative system the oscillations of both kinds of energy 
are sinusoidal. They occur 180° out of phase with each other with 
double the frequency wo of natural oscillations. Both kinetic and po- 
tential energies oscillate between the minimal zero value and the same 
maximal value Frot about the mean value, which equals one half of 
the total energy. All this can easily be seen if you choose the item 
“Energy transformations” in the simulation program (for the case of 
the absence of friction). 


1.3.1.7 The Phase Trajectories with the Same Energy. Consider the os- 
cillations of a conservative oscillator at different initial conditions but with the 
same total energy. What differences do you observe in the plots and the phase 
trajectories in these cases? 


Oscillations with the same total energy in a conservative system are 
characterized by the same phase trajectory. The motion of the repre- 
sentative point in the phase plane and the graphs of such oscillations 
differ only by the time origin. 
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1.3.2 Damped Free Oscillations 


1.3.2.1 The Sequence of Maximal Deflections. Under the action of a weak 
force of viscous friction, the sequence of maximal deflections of a free, damped 
linear oscillator forms a decreasing geometric progression: each consecutive max- 
imal deflection is smaller than the preceding one by the same factor, exp(—yTo) ~ 
1- yTo. 

(a) Calculate the value of the quality factor Q at which the amplitude halves 
during every two complete oscillations. 


The number NV; /2 of full cycles during which the amplitude is halved 
is (In2/7)Q. To obtain Nj/2 = 2, the quality factor Q must be 
27 /In2 = 9.065. 


(b) Input this value in a computer experiment and verify the theoretically pre- 
dicted constant ratio of successive maximal deflections. Note that this ratio does 
not depend on the initial conditions. 

(c) Evaluate the increment of the period of oscillations at this value of the 
quality factor with respect to the period To in the absence of friction (in percent). 
Can you detect the increment in the simulation experiment? The marks on the 
time axis correspond to integer numbers of periods Ty = 27 /wọo without friction. 


To evaluate the period of damped oscillations T = 27/w , we can 
use the following expression: 


| 2 | 1 1 


(1.15) 
For the relative increment of the period AT /To = (T — To)/To we 
find 1/(8Q?) = 0.0015, or 0.15%. 


1.3.2.2* Maximal Deflection after an Initial Push. Imagine, that we excite 
oscillations with an initial push which imparts an initial angular velocity of 2wo 
to the flywheel in its equilibrium position. 

(a) Calculate the first maximal deflection of the flywheel after the excitation 
for the quality factor Q = 5. 

(b) What will be the value of the subsequent extreme deflection which occurs 
in the direction opposite to the first? Verify your answers using the simulation. 


(a,b) For excitation of the oscillator by an initial push from the equi- 
librium position, the initial deflection is zero (y(0) = 0), so that in 
the general solution y(t) = Age~” cos(w t + ôo) we can assume 
ôo = —1/2. Thus for the time dependence of y we find: 


p(t) = Ape sinu,t & Ape ™ sin wot. (1.16) 
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Here we can neglect the difference between w and wo, because at 
Q = 5 the relative decrement of the frequency caused by friction is 
only 0.5%: 


Aw/wo = V1 — 1/(2Q)? — 1 ~ 1/(8Q?) = 0.005. (1.17) 


Assuming w1 = wo, we allow an absolute error in the value of phase, 
wt, which is less than 0.01 if the time t is about T/4. In general, this 
error in phase is not small (i.e., Awt ~ 1) only for those time instants 
t, for which the amplitude already almost damps out: exp(—7yt) ~ 
exp(—y/Aw) © exp(—4Q) « 1. 


Differentiating Eq. (1.16) with respect to time yields the angular ve- 
locity: 
p(t) = Ape (wo cos wot — y sin wot). (1.18) 


To determine the value of Ag in Eq. (1.16), we must use the sec- 
ond initial condition (0) = 2wo. Substituting this condition into 
Eq. (1.18), we find Ap = Q/wo = 2. 


To calculate the first and the second maximal deflections of the fly- 
wheel, we must determine the instants tı and t2 when the angular ve- 
locity becomes zero. From Eq. (1.18) we see that tan wot, = wo/7. 
Since for relatively weak friction wo is much greater than y (in our 
case wo/y = 2Q = 10), the value wot, is close to 7/2. In other 
words, we can assume that the first maximal deflection of the fly- 
wheel occurs at t4 = 7/(2wo), when sin wot in Eq. (1.16) equals +1, 
and the second maximal deflection (to the opposite side) occurs at 
t2 = 37/(2wo), when sin wot in Eq. (1.16) equals —1. (These values 
are exact for zero friction.) Substituting Ag = 2 and these values of 
t into Eq. (1.16), we find the corresponding maximal deflections y1 
and y2: 


pı = 2exp(—y7/2wo) = 2exp(—7/4Q) = 1.71 rad = 98°; 
p2 = —2exp(—3ym/2wo) = —2 exp(—37/4Q) = 
— 1.25 rad = —71.5°. 


1.3.2.3** Complex Initial Conditions. 

(a) Let the initial deflection of the torsion pendulum be 155 degrees, and the 
initial angular velocity be 2w 9. The quality factor Q = 5. Calculate the maximal 
deflection of the flywheel. 

(b) With the same initial deflection (155 degrees) and the same quality fac- 
tor Q = 5 as in the preceding item (a), calculate the maximal deflection of the 
flywheel, if the initial angular velocity equals —2wo. 
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(a, b) When the initial deflection and the initial velocity both have 
non-zero values, it is convenient to use the following form of the 
general solution: 


y(t) = e7 (C coswot + S sin wot). (1.19) 


Just as in the previous problem 1.3.2.2, we neglect here the difference 
between w1 and wo, because for Q = 5 the relative decrement of the 
frequency caused by friction is only 0.5%. 


Differentiating Eq. (1.18) with respect to time yields the following 
expression for the angular velocity: 


(p(t) = eT [(woS' — yC) coswot — (wo + 7S) sinwot]. (1.20) 


The values of constants C and S are determined from the initial con- 
ditions y(0) = yo, (0) = Q: 


C= po, S = (Q + ypo) /wo. 


The instant tı of maximal deflection from the equilibrium position is 
determined by the condition ġ(t1) = 0. Using Eq. (1.20), we find: 


Q 1 
wopo +7Q/wWo  powo/N + 1/(2Q) 


Substituting the values Q = 5, yo = 155° = 2.7 rad, and Q = 2wọ, 
we find that wot; = 0.60. From Eq. (1.19) using the appropriate 
values for C and S, we find the maximal deflection Ymax = 3.31 rad 
= 195°. For the opposite direction of the initial velocity Q = —2wọ 
(part b), we find wot; = 2.47, and Ymax = —2.5 rad = — 143°. 


tan wotı = 


(c) Let the initial deflection of the torsion pendulum be —155 degrees. What 
initial angular velocity would ensure the maximal deflection of 155 degrees (to 
the opposite side), if the quality factor Q = 20? 


We can use considerations based on the conservation of energy to 
evaluate the required value of the initial velocity. 


In this case the magnitude of the maximal deflection in the opposite 
direction must be the same as the initial deflection (155°). So the 
potential energy of the spring at the instant of maximal deflection is 
the same as at the initial moment, and the initial kinetic energy must 
be equal to the energy dissipated during the motion of the flywheel. 


The motion from the initial moment up to the maximal deflection to 
the opposite side lasts approximately one half of the natural period: 
t = Ty /2. We can evaluate the energy dissipated during this motion: 


AE/Ey =1—e7?% =1 -e =1-e°7/@ x n/Q. (1.21) 
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It follows from Eq. (1.21) that for Q = 20 approximately 16% of 
the initial energy Eo is dissipated by friction. Evaluating AF, we 
can neglect the contribution of kinetic energy in the initial energy and 
consider Ey ~ 4 Do. Thus, AE ~ Eor/Q = ¿Dogr/Q. Equat- 
ing this dissipated energy AF to the initial kinetic energy iJ Q?, we 


find: 
T 
wo = 4/548 ~ 0A lol. 


For |pọ| = 155° = 2.7 rad the initial velocity Q must be approxi- 
mately 1.1 wo. 


1.3.2.4* The Phase Trajectory of Damped Oscillations. The phase trajec- 
tory of damped free oscillations for Q > 0.5 is a spiral which makes an infinite 
number of gradually shrinking loops around the focus located at the origin of the 
phase plane. This focus corresponds to the state of rest in the equilibrium position, 
and the phase trajectory approaches it asymptotically. 

(a) How does the size of these loops change while the curve approaches the 
focus? 


The shape of the loops of the phase trajectory does not change as the 
curve approaches the focus (figure 1.1,a). Each subsequent loop is 
similar to the preceding one. All the dimensions (including the mean 
radius) of a loop become smaller with each turn by the same factor 
exp(—qTo) = exp(t/Q) ~ 1 — T/Q (the latter expression is valid 
for Q > 1). 


Figure 1.1: Phase portrait of an underdamped (a, y = 0.1wo) and of an over- 
damped (b, y = 1.1w) linear oscillator with viscous friction. 


(b) Does the time interval during which the representative point makes one 
revolution along the spiral change as the loops of the curve shrink? 
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The time interval during which the representative point makes one 
revolution of the spiral does not change as the loops of the curve 
shrink. This interval equals the period T = 27/w , ~ 27 /wo of the 
periodic (sinusoidal) function in the general solution which describes 
damped free oscillations. 


1.3.2.5* The Dissipation of Energy. Compare the transformation of potential 
energy into kinetic energy (and vice versa) for free undamped oscillations in the 
absence of friction with that for free damped oscillations in the presence of viscous 
friction. 

(a) Show, using a simulation experiment, that if Q = 18.1, the amplitude is 
halved during four complete oscillations, and the total energy is halved during 
two complete oscillations. 

(b) Why is the dissipation of mechanical energy nonuniform during one cycle 
of oscillations? At what instants during a cycle is the time-rate of energy dissipa- 
tion greatest and at what instants is it smallest? 


The instantaneous rate of energy dissipation dE /dt equals the prod- 
uct of the frictional torque (which is proportional to angular velocity 
of the flywheel ) and the angular velocity. Hence dE/dt is pro- 
portional to the square of the angular velocity and, consequently, to 
the momentary value of the kinetic energy. The tangent to the time- 
dependent graph of the total energy is horizontal (instantaneous rate 
of energy dissipation dE’ /dt equals zero) at the moments of maximal 
deflection, at which the flywheel changes its direction of rotation. 


1.3.3 Non-oscillatory Motion of the System 


When viscous friction is strong (Q < 0.5), a disturbed system returns to the 
equilibrium position without oscillating. In the computer simulation, the needle 
asymptotically approaches the zero point from one side. 


1.3.3.1* Non-oscillatory Motion at Critical Damping. Consider the case of 
critical damping, y = wo. 

(a) Why is critical damping preferable in measuring instruments using a needle 
as an indicator? How might your answer apply to the suspension system in an 
automobile? Show that the value of Q in the case of critical damping is 0.5. 


Damping in the suspension system of an automobile prevents the an- 
noying oscillation of the chassis after the car hits a bump or a pothole 
in the road. In measuring instruments such as galvanometers, damp- 
ing is introduced intentionally in order to overcome the problem of 
taking a reading from an oscillating needle. A critically damped sys- 
tem comes to rest in the equilibrium position after an initial distur- 
bance more quickly than in any other case (i.e., than at any other value 
of the damping constant + for the same value of wo). From Eq. (1.8) 
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we find that the value of the quality factor which corresponds to the 
critical damping (y = wo) is Q = 0.5. 


(b) Calculate the maximal angle of deflection if the system, with Q = 0.5, 
receives an initial velocity Q = 5wo in the equilibrium position. 

(c) In what lapse of time does the needle move towards this extreme point? 

Verify your answers by simulating the experiment on the computer. Note that 
the needle approaches the equilibrium position from one side—it does not cross 
the zero point of the dial. 


(b, c) For excitation from the equilibrium position by an initial shove 
which imparts the initial angular velocity Q to the flywheel of the 
critically damped system, the dependence on time of the angle of 
deflection and of the angular velocity in the subsequent motion is 
described by the following functions: 


olt) =Ore™, p(t) =N(1—-yt)e™. (1.22) 


At an instant tı of maximal deflection from the equilibrium position, 
the angular velocity (t,) is zero. From Eq. (1.22) we see that ù = 
0 at yt = 1, that is, at tı 1/y = 1/wo = To/27. Thus for 
critical damping, the flywheel, after the initial shove, moves from the 
equilibrium position to the point of maximal deflection during the 
time interval Ty)/27 = 0.1675, where To is the period of natural 
oscillations in the system in the absence of friction. Substituting this 
value of t in y(t) given by Eq. (1.22), we find the maximal deflection 
from the equilibrium position: Ymax = y(ti) = (Q/7)e7!. For 
Q = 5wo this yields Ymax = 5/e = 1.839 rad = 105.4°. 


1.3.3.2 Critical Damping. 

(a) Prove that the value Q = 0.5 (y = wo) is really critical. Do so by showing 
that at slightly greater values of Q, the needle of a perturbed oscillator executes 
heavily damped oscillations, slowly moving to and fro across the zero point of the 
dial. 

(b) For a critically damped system, express the constants C and C% in the 
general solution y(t) = (Cit + C2) exp(—7t) to the differential equation (for the 
case of critical damping) in terms of the initial displacement (0) = yo and the 
initial angular velocity (0) = Qo. 

For a critically damped system, the constants Cı and C% in the gen- 
eral solution y(t) = (Cit + C2)e~™ to the differential equation are 
expressed in terms of the initial displacement y(0) = (yo and the 


initial angular velocity (0) = © as follows: 


Cy =Q + ypo, C2 = po. (1.23) 


https://gioumeh.com/product/simulations-of-oscillatory-systems-solution/ 


1.3. PROBLEMS WITH ANSWERS AND SOLUTIONS 17 


Specifically, when motion of the system is excited from the equilib- 
rium position (pọ = 0) by a push which imparts to the flywheel an 
initial angular velocity Q, we get C1 = Q, Cz = 0. The subsequent 
motion is described by the function y(t) = Q texp(—yt). 


(c) Is it possible for a critically damped system to move after an initial distur- 
bance according to a pure exponential law? If so, what initial conditions give rise 
to such motion? What is the phase trajectory of this motion? Prove your answers 
experimentally. 


The motion of the flywheel after excitation is described by a pure ex- 
ponential function of time if in the general solution C4 is zero. From 
Eq. (1.23) we see that this situation occurs if the flywheel, being dis- 
placed from the equilibrium position through an arbitrary angle yo, is 


given the initial angular velocity Q = —yo. For such initial condi- 
tions both the angle y and the angular velocity p decrease exponen- 
tially: 


y(t) = po exp(—t), p(t) = Qexp(—7t). 


The angular velocity ~ at any instant is proportional to the angle of 
deflection: ¢ = —v~yy. In the phase plane, this equation gives a 
straight line which passes through the origin with a negative slope 
—y. The representative point moves along this line from the initial 
point (Yo, 22), approaching the origin asymptotically. The duration of 
the motion is characterized by an interval r = 1/7. 


(d) At what initial conditions the flywheel of a disturbed critically damped 
system will cross the equilibrium position? For a given initial displacement yo, 
what initial angular velocity Q should you impart to the flywheel of the critically 
damped oscillator in order it crossed the equilibrium position after a lapse of time 
t = 37, where To = 2r /wo is the natural period (the period of oscillations in the 
absence of friction)? 


The flywheel of a critically damped system crosses the equilibrium 
position after a disturbance if initially it is displaced from the equi- 
librium position and receives a strong enough initial shove directed 
towards the equilibrium position. For a given initial displacement yo 
the initial angular velocity Q must have a sign opposite that of the ini- 
tial displacement and its magnitude must exceed the value |yyo| = 
[wopo]. 

In the phase plane any phase trajectory that starts above the straight 
line p = —yọ approaches the node (the origin of the phase plane) 
moving along this line from the side of positive values of y. There- 
fore, for a negative initial displacement pọ < 0, in order to cross the 
equilibrium position y = 0 (the ordinate axis of the phase plane), the 
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phase trajectory must start from a point lying above the straight line 
py = — yy. Thus, Q > —yyo. On the other side, for yo > 0 the ini- 


tial point of the phase curve must lie below the line p = — yọ. That 
is, Q < —yo. For example, if the initial deflection of the flywheel is 
—20° = —0.35 rad, we must set a positive value of the initial velocity 


which is greater than 0.35 y = 0.35 wo. 


If the flywheel is initially deflected through an angle Yo and is given 
the critical initial velocity Q = —yyo = —wWo Po, its motion towards 
the equilibrium position lasts indefinitely. To reach the equilibrium 
position after a finite lapse of time, the initial velocity must be greater 
in magnitude than the critical value: Q = —woyo(1 + 6). We can 
find the fractional surplus 6 of the initial velocity that is necessary 
for reaching the equilibrium position during a given time (say, during 
t = 379), using Eq. (1.23) for constants C and Cy. In this case 
C2 = go, Ci = —woyod, and the angle of deflection depends on 
time as follows: 


p(t) = poll — dwot)e*™. 


Thus, the instant ¢ at which the flywheel crosses the equilibrium po- 
sition equals 1/(dwo) = To/(275). The lapse of time is 375 for 
ð = 1/6 ~ 0.053. The required value of ô is independent of the ini- 
tial deflection. For a given initial displacement, say —20° = —0.35 
rad, we should enter an initial velocity whose magnitude is 5.3% 
greater than the critical value 0.35wo: Q œ~ 0.37wo. Then the fly- 
wheel crosses the equilibrium position at t ~ 37. 


1.3.3.3* Motion of an Overdamped System. 

(a) For arbitrary initial conditions (y(0) = Yo, ¢(0) = Qo), express the values 
of Cı and C% in the general solution to the differential equation for an overdamped 
system in terms of Yo and Qo. 


The general solution to the differential equation of the overdamped 
system is a superposition of two exponential functions: 


y(t) = C1e™t + Coe", where a12 = —y+ 4/7? — w2. (1.24) 


Since both a; and az are negative, these exponential functions de- 
crease in magnitude with time t. Constants C and C% in Eq. (1.24) 
are determined by the initial conditions y(0) = yo, g(0) = Q: 


Q — azYyo aiıpo- NQ 


Cy = = 


oa y= we’ 2 2 _ 
0 Y wo 


(1.25) 


(b) At what initial conditions the motion of an overdamped system will be 
described by a monoexponential function of time? What are the phase trajectories 
that correspond to such motions? 
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Motion of the disturbed system is described by a single exponential 
function if either C; = 0 or C2 = 0. In the first case the asymptotic 
exponential approach to the equilibrium position y(t) = woe? is 
characterized by a time constant 72 = —1/a2. This case occurs if 
the initial conditions satisfy the requirement Q = a2gy~o. The second 
case y(t) = ye" is characterized by a smaller time constant 7, = 
—1/a, and occurs if the initial conditions satisfy the requirement 
Q = aro. 


The phase trajectories for these cases are straight lines passing through 
the origin of the phase plane and having negative slopes (see fig- 
ure 1.1), b). Their equations are ý = ayy and ý = agy respectively. 


(c) Explain, why at arbitrary initial conditions non-oscillatory motion of the 
flywheel towards the equilibrium position occurs more slowly and requires more 
time than at critical damping. Is it possible for an overdamped system to return 
to the equilibrium position faster than for the critically damped system with the 
same wọ? If so, what conditions of excitation ensure the motion? 


For arbitrary initial conditions, when the constant C in Eq. (1.24) 
is nonzero, the duration of the process of exponential approach to 
the equilibrium position is characterized by a time 7, = —1/a, = 
1/(y — y7? — we). This time is greater than the characteristic time 
rT = 1/7 of a critically damped system. The only exception is the 
case for which C = 0. In this case the duration of the exponential 
motion towards the equilibrium position is characterized by a time 
T2 = —1/az = 1/(y + y7? — wg), a duration that is shorter than 
T = 1/4. This case of relatively fast approach of the system to the 
equilibrium position is realized only under specific initial conditions, 
namely for Q = azYvo. 


(d) What is the principal difference between the phase trajectories correspond- 
ing to a non-oscillatory motion and those corresponding to damped oscillations? 


The principal difference between phase trajectories that correspond 
to non-oscillatory motion and to damped oscillations is in their topol- 
ogy. Phase trajectories of damped oscillations approach the origin of 
the phase plane (the state of rest in the equilibrium position) winding 
around it. The spiral makes an infinite number of turns around the 
focal point. Phase trajectories of an overdamped system approach the 
origin of the phase plane directly, without spiraling. At this nodal 
point the curves have a common tangent ý = a y. The tangent is 
one of these phase trajectories. Motion along it occurs for initial con- 
ditions of the kind Q = a iyo. The only phase trajectory for which 
this line is not a tangent is the straight line p = agy. It corresponds 
to the fastest damping of exponential motion, as discussed above. 
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(e) Is it possible for an overdamped system (y > wo) to cross the equilibrium 
position after excitation? If so, what initial conditions give rise to such motion? 
Is it possible for the oscillator to cross the equilibrium position more than once? 


A critically damped or an overdamped system can cross the equilib- 
rium position after an initial excitation only once. This situation oc- 
curs if initially the flywheel is displaced from the equilibrium position 
and is given a large enough initial angular velocity directed towards 
the equilibrium position. From the phase portrait of an overdamped 
system we can conclude that all the phase trajectories which start at 
a point lying above the more inclined straight line 2 (g = agy) ap- 
proach the node from the side of positive values of y, and vice versa. 
The system in motion crosses the equilibrium position if the phase 
trajectory intersects the ordinate axis y = 0 of the phase plane. Con- 
sequently, if the flywheel is initially deflected to the left (pọ < 0), in 
order to cross the equilibrium position, its initial state must be repre- 
sented by a point lying above the line 2. Thus, at yo < O the initial 
angular velocity Q must be positive, and its magnitude must exceed 
the value a290 = (y + y7? — w) |yol. If the flywheel is deflected 
to the right (po > 0), the initial velocity must be negative (directed 
towards the equilibrium position), and its magnitude must satisfy the 
same condition as in the previous case. 
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